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Semiclassical gravity is investigated in a large class of asymptotically flat, static, spherically 
symmetric spacetimes including those containing static stars, black holes, and wormholes. 
Specifically the stress-energy tensors of massless free spin and spin | fields are computed to 
leading order in the asymptotic regions of these spacetimes. This is done for spin fields in 
Schwarzschild spacetime using a WKB approximation. It is done numerically for the spin ^ 
field in Schwarzschild, extreme Reissner-Nordstrom, and various wormhole spacetimes. And 
it is done by finding analytic solutions to the leading order mode equations in a large class 
of asymptotically flat static spherically symmetric spacetimes. Agreement is shown between 
these various computational methods. R is found that for all of the spacetimes considered, 
the energy density and pressure in the asymptotic region are proportional to r -5 to leading 
order. Furthermore, for the spin 1/2 field and the conformally coupled scalar field, the 
stress-energy tensor depends only on the leading order geometry in the far field limit. This 
is also true for the minimally coupled scalar field for spacetimes containing either a static 
star or a black hole, but not for spacetimes containing a wormhole. 
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I. INTRODUCTION 

A number of studies have been carried out of the leading order quantum corrections to the 
gravitational potential due to quantized fields in static spherically symmetric spacetimes in the 
weak field limit [liul]. In two cases the leading order behaviors of the energy densities 

and pressures of the quantized fields were also computed. In this limit it is assumed that the 
gravitational field is everywhere weak such as in the case of a planet or a nonrelativistic star. The 
far field limit in an asymptotically flat spacetime is similar to the weak field limit in that, in the 
asymptotic region, the spacetime is nearly flat. The difference is that the spacetime curvature 
need not be small everywhere. Since quantum field theory is intrinsically nonlocal, it is in principle 
more difficult to study quantum effects in the far field limit when the spacetime curvature is not 
small everywhere and the fields are in the vacuum state which is natural for the asymptotically flat 
region, often called the Boulware state [8]. 1 In fact there are only two previous calculations that we 
are aware of which focus on the far field limit: Anderson, Balbinot, and Fabbri [lo| computed the 
leading order asymptotic behavior for the stress-energy tensor, (T a j,) for massless scalar fields with 
arbitrary coupling £ to the scalar curvature in the Boulware state in Schwarzschild spacetime; and 
Garbarz, Giribet, and Mazzitelli calculated for the same fields in D dimensional Schwarzschild- 
Tangherlini spacetime with D > 4 the leading order asymptotic behaviors of both ((p 2 ) and (T a h). 

An interesting question that can be addressed by comparisons of calculations in the weak field 
and far field limits is the extent to which leading order quantum effects are local when the fields 
are in the Boulware vacuum state. By "local" we mean that in the far field limit the leading order 
terms in the stress-energy tensor depend only on the leading order deviations of the metric from 
a flat space metric. Comparison of the calculations in Ref. 10] for Schwarzschild spacetime with 
the weak field calculations of Ref. [J] for nonrelativistic static spherically symmetric stars shows 
that while there is agreement in the leading order asymptotic behaviors of the energy density and 
pressure for massless scalar fields with conformal and minimal coupling to the scalar curvature, for 
all other couplings there are differences. Further there are also disagreements in the leading order 
behavior of the quantity (<p 2 ) for all couplings to the scalar curvature except minimal coupling. In 
Ref. pjj two of us computed the leading order difference in the stress-energy tensors in the far field 
limit for a Schwarzschild black hole and a static spherically symmetric star for massless scalar fields 



1 This statement does not apply to thermal states such as the Hartle-Hawking-Israel state [{J, because for such states 
the leading order behavior of the stress-energy tensor for a quantized field in the far field region is the same as it 
would be in Minkowski space. 
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in the Boulware state. It was shown that in most cases the leading order behaviors of the quantities 
(<p 2 ) and (T a b) have contributions from the local geometry as well as nonlocal contributions from 
a zero frequency and zero angular momentum mode. However, when the curvature coupling is 
minimal there is no difference at leading order between the case when a black hole is present and 
the case when a star is present. When the coupling is conformal there is a difference for (</> 2 ), but 
for (T a b) the differences cancel, presumably because of the conformal symmetry. In Ref. [5j the 
calculation in Ref. was generalized to the case of D dimensions with D > 4. As for D = 4, it 
was found that there is no leading order difference in (eft 2 ) and (Tab) f° r minimal coupling and no 
leading order difference in (T a b) for conformal coupling in D dimensions. 

In this paper, analytic solutions to the leading order mode equations in the far field limit of 
a large class of asymptotically flat, static, spherically symmetric spacetimes are found for both 
massless spin and spin 1/2 fields. From these solutions analytic expressions for (cf 2 ) for scalar 
fields and (T a b) for both types of fields are derived when the fields are in the Boulware vacuum 
state. For the spin 1/2 field it is shown that the leading order behavior of (T a b) depends only on 
the geometry in the far field region. For the conformally coupled scalar field it is shown that the 
leading order behavior of (T a b) also depends only on the geometry in the far field region, but this 
is not true for (4> 2 ). For the minimally coupled scalar field it is shown that there are cases, such 
as spacetimes containing a star or a black hole, where the leading order behaviors of both (cj) 2 ) 
and (T a b) depend only on the geometry in the far field region. However, for the minimally coupled 
scalar field in wormhole spacetimes the leading order behaviors of these quantities depend upon 
the entire geometry. 

Numerical computations of the full stress-energy tensor for the spin 1/2 field in Schwarzschild, 
extreme Reissner-Nordstrom, and three wormhole spacetimes are analyzed and displayed. These 
provide important checks on the analytic calculations. In each case the leading order behavior in 
the far field limit has been found to agree with the analytic results to within the numerical error 
of the computations. 

In Section |H] the computation of the stress-energy tensor for a massless scalar field is reviewed. 
In Section IIIIl a somewhat peculiar method mentioned in [l(| of using the WKB approximation 
and conservation of the stress-energy tensor to compute the leading order behavior of the stress- 
energy tensor for massless scalar fields in the far field limit of Schwarzschild spacetime is given. 
Section IIVI contains a review of the computation of the stress-energy tensor for the massless spin 
1/2 field. In Section [V] an analytic computation is made of the leading order behavior of the stress- 
energy tensor for the massless spin 1/2 field in the far field limit. A similar analytic computation is 
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made for a massless scalar field with arbitrary coupling £ to the scalar curvature in Section IVII In 
Section fVIII some of the numerical results for the massless spin 1/2 field are displayed and discussed. 
Section IVIIII contains a comparison between our results and previous analytic calculations of the 
leading order quantum corrections to the gravitational potential in Schwarzschild spacetime. Some 
conclusions are given in Section IIXI Various identities relating to modified Bessel functions which 
were used in the analytic computations are proven in the Appendix. Throughout we use units such 
that h = c = G = 1, and our conventions are those of Misner, Thorne, and Wheeler [l^ . 



II. GENERAL FORM FOR THE STRESS-ENERGY TENSOR FOR MASSLESS 

SCALAR FIELDS 

In this section the method developed to compute the stress-energy tensor for scalar fields in 
static spherically symmetric spacetimes given in Ref. [jjj is adapted to the special case of computing 
the leading order components of this tensor for a class of asymptotically flat spacetimes in the far 
field limit when the fields are in the Boulware state. 

The metric for a general static spherically symmetric spacetime can be written as 

ds 2 = -f{r)dt 2 + h{r)dr 2 + rW . (2.1) 

In this paper we consider the class of asymptotically flat spacetimes for which the metric functions 
in the large r limit have the behaviors 

f(r) = l- — + o(\) and h(r) = 1 + — + O ( \ ), (2.2) 

where A and B are two parameters describing the metric at large distances; for example, for both 
the Schwarzschild and Reissner-Nordstrom metrics, A = B = M. 

* Re, Q it wa S show,, tt* «,e s tr e„ t e„ SOr for a sca laI fl e ld in a g e„e r a! s t a ti c 
spherically symmetric spacetime can be written in terms of an analytic term plus a term involving 
sums and integrals over five different combinations of the modes, which must often be compute 
numerically, 

{Tab) = {Tab) numeric + {Tab) analytic • 

(2.3) 

For fields in the Boulware state, the quantity (T a f ) ) ana i yt i c behaves as r -6 at large r, but as we shall 
see, (T a 6) numer i c dominates it and behaves as r -5 . For a massless scalar field with coupling £ to 
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the scalar curvature, (T a f ) )numeric has the asymptotic form 2 



<T^_ ri c = (2£ + l)S 1+ (2£ - |) S 2 + (2£ - |) - (2£ - i)^5 5 



(T 6 



r /numeric 



1 2£ 1 

~Si + — TT^Sz H #4 + Tr - ?^ > 



2r 2 



8r 2 



(2^-i)5 1 + (2C-i)5 2 + ^3-^4-^5 



(2.4a) 
(2.4b) 
(2.4c) 



Here the 5 n 's are sums and integrals over the radial modes p w (,[r) and q u &{r) of the Euclidean 
Green function. In the Boulware vacuum state we anticipate that to leading order the nonzero 
components of the stress-energy tensor in the asymptotically flat region go like 0(r -5 ). Hence we 
need to determine S\ and S2 to 0(r -5 ), S3 and £5 to 0(r -3 ), and S4 to 0(r~ 4 ). To this order we 
find 3 



dhJ 

4^2 



E 



.e=o 



{2i+l) Pujl {r)q ut {r) - \ - -3 



+ w + 



2Aj 



E 

.£=0 



du> 
4^2 



4vr 2 ] 

.£=0 



ar ur \ 2r 2r^ 



1 A + 2B\ B-A 
; 1 + 2r 4 



,,3 , , , , (^+|) 2 (^+^) 2 A co 2 r 3Alu 2 



2r 2 w 3 8vlra; 3 a; Aw _Bw 
3 3 4 6r 3r 



5, 



«=0 



4^ 2 

/•°° dj 
4^ 



(2^+1)^^(0^(0] + ^ + ^ 
(21 + D/o ('•)'/.,/! '•) - - - 4 + + 



2Aj 
2Aj] 



The same type of split also occurs for (<j> 2 ) with 



(2.5a) 



AA + 2B\ uj 3 A + B , 



(2.5c) 
(2.5d) 
(2.5e) 



/ numeric 



(2.6) 



For fields in the Boulware state, (0 2 ) an aiytic goes like r 4 at large r so, as for the stress-energy 
tensor, ((/> 2 ) n umeric dominates. 



2 The complete forms for (T a f,) n umcric and the S n sums are given in Ref. [lj]. Here we present them only to leading 
order in the metric perturbations A and B. 

3 Note that the dependence on the mode functions p^e and q u t in these expressions is exact. It is the subtraction 
terms which are approximate. 
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The quantities Pui(r) and q^eif) are radial mode functions for the Euclidean Green function. 
They satisfy the equation 

d 2 S 



2 1 df 1 dh\ dS 
dr 2 \r 2f dr 2h dr J dr 



co 2 h {£ 2 + e)h 

f r 2 



S = 



(2.7) 



where R is the scalar curvature. To first order in the metric parameters A and B, this equation is 



+ - + 



A + B\ dS 



dr 2 

Also note that the relationships 

Si 
S 2 



dr 



u; 2 [ 1 



2A + 2B 



S 



+ • 



1 + ^)5 = 0. 



dSs 
dr 
1 

2 



2h 



f 



S 



2h 



S 3 + 



h'\ 



2 / 

r + 2f 2hJ 



dS<i 
dr 



2£Rh 



2r 2 



(2.9a) 
(2.9b) 



are exact. 4 

The modes are normalized so that they satisfy the Wronskian condition 



/ sdquiir) dp ui (r) 1 / h 

P^—jz —^{r) = (j 



1/2 



(2.10) 



dr dr 

In an asymptotically flat spacetime one boundary condition is that q u i{r) is finite in the limit 
r — > oo. The boundary condition for p^(.{r) depends on the behavior of the geometry at small 
values of r. If the geometry is regular at r = 0, such as in the case of a static star, then p w i{r) 
is regular there as well. If there is an event horizon then p^e( r ) should be regular at the event 
horizon. For a wormhole p^ii^r) should go to zero in the asymptotically flat region on the other 
side of the wormhole. 



III. WKB APPROXIMATION FOR THE STRESS-ENERGY TENSOR FOR A SCALAR 

FIELD 

As shown in Ref. [13], for a scalar field the WKB approximation for the radial mode functions 
is obtained by first making the transformation 

A* = (2,2^)1/2 ex P {/ W (j) *■} . ( 3 - la ) 

1 f r^fh^ 1 ' 2 A 



= (2^)1/2 eXp W W \j) dr ( ■ (3 - lb) 



4 To see this use the exact form of the mode sums in [13J along with the mode equation (|2.7p . 
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Note that the Wronskian condition (|2.10|) is identically satisfied by these expressions. Substitution 
into the mode equation (12.7[) gives in both cases 



W 1 



n 2 + Vi + v 2 + 



1 



/ d 2 W d //\ 1 dW 3/ / 1 dW\ 



hW dr 2 dr\h) 2W dr 2h \W dr J 



(3.2) 



with 



n 2 

V-i 

v 2 



1 df f dh f 



2rh dr 2rh 2 dr 4r 2 



(3.3a) 

(3.3b) 
(3.3c) 



Eq. (|3.2p can be solved iteratively. The zeroth order solution is W = VI. The second order solution 
is 



W = £l + —(V 1 + V 2 ) 



m 3 1 + 4 



/ d 2 Q d ff^ 1 dfl 3/1 (dfl x 2 



hQ 2 dr 2 



(3.4) 



dr \h J 2tt 2 dr 2 h Q 3 \ dr 

An approximation to (T a b) numer i c can be obtained by first substituting Eqs. (|3.1|) into the sums 
in Eqs. (]2.5p . Approximations for these sums, which we call (SVOwKBfin, can be obtained by substi- 
tuting the WKB expansion to some order for W and its derivatives into the resulting expressions, 
keeping only terms which are smaller than or equal to the given order of the WKB expansion. 
Note that the order is increased by one for every radial derivative of W. The details of how the 
mode sums and integrals are computed are given in Appendix F of Ref. The results are then 
substituted into Eqs. (12. 4p to obtain the approximation which we call (T a ;,)wKBfin- While the exact 
expression for (T a b) numer i c is conserved, in general (T n ())wKBfin is not a conserved tensor. 

If a second order WKB expansion is used for Schwarzschild spacetime, where 



/ = h' 1 = 1 - 2M/r , 



(3.5) 



then at lowest order there are terms in the expressions for (T a (,)wKBfin which at large r go like 
1/r 4 , M/r 5 , M 2 /r 6 , and so forth. 5 The lowest order terms are present in the M — > limit which 
is flat space. Thus there is some justification in ignoring them since the full renormalized stress 
tensor is zero in flat space. The terms of order M 2 /r 6 and higher are of subleading order at large 
r, so it is reasonable to ignore them. That leaves the terms of order M/r 5 . Terms of this form do 
not appear in (T ,a fe) ana i y ti c , which has instead leading order terms proportional to M 2 /r 6 . 



5 Higher order terms in the WKB expansion also generate terms of order 1/r 4 , M/r 5 , etc.; however, the correct 
leading order behaviors for the sums S3, 5*4, and 5*5 are obtained by using only a second order WKB expansion. 
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The terms of order M/r 5 in the stress-energy tensor come from terms of order M/r 3 in 
(5 , 3)wKBfin and (S , 5 )wKBfin, M/r 4 in (5 4 ) W KBfin and M/r 5 in (5i)wKBfin and (S 2 )wKBfin- We 
find that the terms of these orders in (S^wKBfin, (^jwKBfin, and (S^wKBfin are 

M 

5 4 = -g^i > ( 3 - 6b ) 

s " = 2W ■ (3 - 6c) 

These are in agreement with both the numerical computations discussed in flo| and the exact 
analytic computations shown in Section VI. This is not the case for the corresponding terms in 
(Si)wKBfin and (S^wKBfin- However, if Eqs. (]3.6p along with expressions of the form S\ = aM/r 5 
and S2 = bM/r 5 , are substituted into Eqs. ([2.4]) . then the result is conserved to leading order if 
and only if 

» - ■ ^ 

The resulting expressions for the order M/r 5 contributions to S± and S2 are in agreement with 
both our numerical and our exact analytical results. If these results are substituted into Eq. (12.41) 
then one finds that the leading order behavior is 

M 

(T\) = [diag (-1, 2, -3, -3) + 5£diag (2, -2, 3, 3)] , (3.8) 



which is the same as that given in Ref. [101 ] . 

Note that one might expect a similar calculation to work for the massless spin 1/2 field. However, 
because the second order WKB approximation, when used in the way described above, does not 
give the correct leading order behaviors for all of the sums in Eqs. (]2.5p . we do not pursue this 
method further here. 



IV. GENERAL FORM OF THE STRESS-ENERGY TENSOR FOR THE MASSLESS 

SPIN 1/2 FIELD 

In this section the method developed to compute the stress-energy tensor for the massless spin 
i field in static spherically symmetric spacetimes given in Ref. is adapted to the special case of 
computing the leading order components of this tensor for a class of asymptotically flat spacetimes 
in the far field limit. 



f. Q 



In Ref. [14| it was shown that the stress-energy tensor for a massless spin i field in a general 
static spherically symmetric spacetime can be written in the same way as Eq. (|2.3|) with an analytic 
term plus a term involving sums and integrals over the modes. As for the scalar field, (T a ;,) ana iyti c 
gives a leading order contribution at large r of 1/r 6 which is less than the 1/r 5 contribution that 
we expect from (r a fe) numeric . For a metric with the asymptotic form (|2.2|) . (T a ;,) numer i C has the 
asymptotic form 6 



{T r 



numeric 



7T^ 



1 



duo 



1 



+ 



2,4 



1 + — 

r ) 



(4.1a) 



r numeric 



/ numeric 



ZS / du 
K JO 



A 



1 



+ -r ) uA 2 - [ ~n +~T ] u 2 A x + ( 1 + ^ 



2A 



4A 
3r 



uj° + 



A + B 



6r" 



-UJ 



1 

2^2 



doj 



1 



+ 



.4 



w4 2 + TT + 



8A 



A + B 



oj ' 



-OJ 



(4.1b) 
(4.1c) 



3 3r J 6r 3 

The asymptotic forms of the sums A\ and A 2 that we use for the calculations in Section [V] are not 
simply the asymptotic forms of the exact equations for A\ and A 2 given in Ref. [l4]. To eventually 
derive the form of the equations which we use, it is useful to first discuss the exact expressions for 
A\ and A 2 given in Ref. [14]. They are 



A i= E C;+i)J2A-Cj-i)G£j-iG£j-i + 7^ 

j=l/2 1 1 



A, 



oo 

£ 

i=i/2 



7 -2_I 



+ 



2 

r oj 

w 



Here the sum is over half integer values of j, with j 



+ g in the notation of Ref. 



(4.2a) 
(4.2b) 



and G p , _■ are radial mode functions for the Euclidean Green function, as are F q ■ and Gj, The 
boundary conditions for the mode functions are the same as those for the scalar field. The notation 
is also similar, so F q , ■ and G q , ■ are the modes which vanish at spatial infinity, while F p ■ and G p , ■ 
are the modes which are regular at the origin if there is a star, or at the event horizon if there is a 
black hole. They vanish at spatial infinity in the other universe if there is a wormhole. Each pair 
satisfies the coupled set of equations 



1 

h 1 / 2 dr u ' j 



OJ 



J + -2 



M :G 



G, 



u,3 i 



, 1/2 ^J-T r 



(4.3a) 
(4.3b) 



6 The complete form for {T a b} numeric is given in Ref. [IJ]. Here we present it only to leading order in the metric 
perturbations A and B. 
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The radial mode functions also satisfy the Wronskian condition 



(4.4) 



Equations (|4.2p can be recast in a more useful form in which the subscripts are consistently the 
same. As a first step, rewrite each sum as the limit of a finite sum, and then shift j — > j+1 on 
those terms in the sum for which the subscript is j— 1 rather than j. The result is 

J r n J-l 



lim 

J— ¥00 




fl/2 



i=-i/2 



An 



lim 

J— >oo 



lim 

J^oo 



lim 

J— >oo 



.5=1/2 1 J 
J 

E 



U1,J U1,J ( ' 



(4.5a) 



J 

E 

.5=1/2 



(i 2 -|)^<5 - ^ + 



2 



2 

r w 
IT 



i=-i/2 



(J 2 +2J+!)< J ^ iJ 



(4.5b) 



Here the Wronskian condition (|4.4p has been used in the final step to simplify the expression for 

The limiting terms are to be evaluated in the large J limit and for this purpose a WKB-like 
expansion can be used for the mode functions. We first define new mode functions ■ and ■ 



as the ratio 



Z 



G, 



Then Eqs. (14.31) combine to become the single equation 



1 d 



1 - + Bz±I Zl 



dr ^ fl/2 ■ 

Using the quadratic formula one finds the formal expression 



(4.6) 



(4.7) 



U+l)f 1/2 



1 ~\ — n . - /„ —T 



1/2 



(4. 



coh l / 2 dr u<3 

This equation can be solved iteratively. The zeroth order solution is obtained by setting dZ^j/dr = 
on the right hand side. Since there are two solutions, one must correspond to Z q and the other 
to Z p . Because the F q and G q mode functions vanish at infinity, their magnitudes must decrease 
at large r. From Eq. ()4.3bp it can be seen that the magnitude of F q can decrease as r increases 
only if G q and F q have opposite signs. Therefore 



z h = p y G h < ° 



(4.9) 
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at large r. Thus Z q must correspond to the solution with the minus sign, and therefore Z p 
corresponds to the solution with the plus sign. In the large j limit, we then find that 



e.. = (2J+1)/1/2 + 0( i). 



Z q • 



oor 
ujr 



(4.10a) 
(4.10b) 



(2j+l)/i/a + ° (r } • 

Any product of p- and q- modes can be written, with the help of Eq. (14. 4p . in terms of the Z's; 
for example, 

1 



fig ftp _ 



9 > 
w,3> 



p9 /oP 



z q . 



1 > 



Substituting the approximate form (|4.10p gives 



(4.11a) 
(4.11b) 



lim (J+^GIjG^j 



J— >oo 
T'2 



lim (J 2 +2J+f)F« G^, 



J— >oo 



2/V2 
o 

r u 
~ 4/ 



Substituting Eqs. (|4.12p into Eqs. (f4"3j) then yields 

OO r- 



A 2 



3=1/2 

£ 

3=1/2 



/ l/2 



+ 



2/1/2 



r cj 



+ 



2 

r oj 

V 



Eqs. (|4.13|) are completely general. To first order in A and B they become 



Ax = E [(i+l)(<i<i-^, i ^ )i )+r + ^ 

3 = 1/2 



A 2 



+ ±r 2 w + i^lrw . 



3=1/2 



(4.12a) 
(4.12b) 



(4.13a) 
(4.13b) 



(4.14a) 
(4.14b) 



It is this form which we use for our computations in Section [Vj 



V. ANALYTIC COMPUTATION OF THE STRESS-ENERGY TENSOR FOR THE SPIN 

i FIELD 

In this section the asymptotic behavior of the stress-energy tensor for a massless spin 1/2 field 
is computed analytically. To do so we first use Eq. (|2.2p in Eqs. (|4.3p to write in the large r limit 



12 



the approximate mode equations 

dF ud j + \ A , B\ „ ^ , A + £ 



1 + - F w j = w 1 + G w j , (5.1a) 



dr r \ r 

It is useful to change variables. First let 



+ 1 + - G uJ = uj 1 + F UJ . (5.1b) 



i + 2 



F U j = F uJ + J -^BF^ 3 + uBG uJ , (5.2a) 

G u j = G u>j - J -^BG u>j + uBFuj , (5.2b) 

and define 

x = ijjr . (5-3) 



The mode equations are then 



dF ud j + f, , Cuj 



F u j = 1 + G UJ , (5.4a) 



dx x ' \ x 

+ J -^G^ = + — (5.1b) 



— j 1 "wj — IH ^ Id J ) 

ax x \ x J 

with 

+ (5.5) 

The functions and Gr w ,- satisfy, to first order in jB/r, the same Wronskian condition (|4.4j) 
that F^j and G^j do. 

If C = 0, it is not hard to show that the solutions to Eqs. (|5.4I) will be given by (x/uj) 1 ^ 2 Ij(x) 
and (x/uj)^/ 2 Kj(x) for and (x/w) 1 / 2 / 7+ i(x) and (x/w) 1//2 -fTj + i(x) for G u j, where Ij and 
are modified Bessel functions. Many of the properties of these functions, including a number of 
useful identities, can be found in Appendix „ The Ij's are well behaved at small x, but diverge at 
large x, while the Kj's are infinite at x = 0, but damp to zero in the limit x — > oo. Recalling that 
x = wr, we thus anticipate that the F q and G q modes will primarily involve the Kj's. 

For C 7^ we can continue to write F U j and G^j in terms of (x/uj) 1 ^ 2 Ij(x) and (x/uj) 1 ^ 2 Kj(x) 
by defining the new pairs of variables aF u P^ j( x ) an d ofl^ -{x), -(x) such that for each pair 

1 /2 

Fujj(x) = (^j [a U) j(x)Ij(x) - P U) j(x)Kj(x)] , (5.6a) 

1 /2 

Gujix) = (^J [a U j(x)I j+ i(x) + P U j(x)Kj +1 (x)] . (5.6b) 
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Using the Wronskian condition (|A.5P for the Bessel functions gives to first order in jB/r the 
condition 



(5.7) 



Substituting Eqs. (]5.6h into the mode equations (I5.4j) and using the recursion relations (IA.3j) to 
eliminate the derivatives acting on the modified Bessel functions, yields coupled linear expressions 
for a' u Ax) and -(x). These can then be simplified using Eq. (lA~5jl to give 

q^.(x) = Cwfo^tofo+iOr)^ , (5.8a) 

fc d (x) = Coj {a^(x) [l]{x) - l] +1 (x)} - p uJ (x) [I J+1 (x)K J+1 (x) + I^K^x)]} , (5.8b) 

where primes denote derivatives with respect to x. 

To compute the stress-energy tensor, we shall require the products of the p- and g-solutions, 
and for this reason it is helpful to work with the following equations for the relevant products of 
the a's and /3's: 



a i,j( x ) a l,j( x ) 



Cu>{2c% tj (x)al d (x) [I j+1 (x)K j+1 (x) +I j {x)K j {x)} 

+ [^MjW+^^W] [K] +l (x)-K](x)]} , (5.9a) 
Coj{plj(x)P q ^(x) [K] +1 (x)-K](x)] 

C U {$>j{x)%j{x) [K] +l (x)-K](x)] 

+ [I](x)-I] +I (x)]} , 

Coj{ [aZ, d (x)fi d {x) + ^(x^x)] [l](x) - l] +1 (x)] 



(5.9b) 



(5.9c) 



(5.9d) 



We want to solve these equations to first order in Coj. The solutions will be substituted into 
Eqs. (I4.14p and the sums over j will be computed. Then the results will be substituted into 
Eqs. (14. ip and the integrals over oj will be computed. For large values of j the condition (15.71) 
breaks down and for large enough values of oj the quantity Coj is not small. Therefore one might be 
concerned that the leading order behavior of the stress-energy tensor will not be correct. However, 
to leading order the nonzero components of the stress-energy tensor go like D/r^ for some constant 
D. If we temporarily consider units in which H = c = 1 but G / 1, then both C and D have 
units of length. Thus D can be proportional to C but it cannot be proportional to C 2 or any 
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other power of C. For this reason it should be sufficient to only keep terms to order C in both the 
condition (I5.7P and the solutions of Eqs. (15. 9|) . 

In flat space the only nonzero functions are -(ac) and -(x), and these are constants. For 
the class of spacetimes we are considering we expect that this will continue to be true to leading 
order in a perturbative expansion. However, we cannot rule out the possibility that pathological 
metrics exist for which this condition would be violated. Thus our calculations may not cover every 
possible spacetime with an asymptotic metric of the form (|2.2p . but as can be seen in Sec. IVIIl 
they definitely apply to two standard black hole metrics and some wormhole metrics and there is 
every reason to believe that they apply more generally to most or all static spherically symmetric 
black holes, static stars, and wormholes. 

Hence, to leading order, we assume that on the right hand side of Eq. (|5.9p all products 
other than a£, j(x)f3^ j(x) vanish. Furthermore, the constraint (|5.7j) demands that to this order 



1, so Eqs. (|5.9p simplify to 

<;(*)<^)]' = Cu[K] +l {x)-K]{x) 



Clo[I]{x)-I 2 0+1 (x)] . 



(5.10a) 
(5.10b) 
(5.10c) 
(5.10d) 



We now need to find suitable boundary conditions for Eqs. (|5.10p . This can be accomplished in 
part by looking at the solutions to Eqs. (|5.8p in the limit x — > oo. Using the usual asymptotic 
expansions for the Bessel functions one can show that for arbitrarily large values of x 



I j+ i(x)K j+1 (x) +Ij(x)Kj(x) 
K] +1 {x)-K]{x) 



-n-l 



I 2 Ax) 



(x) 



n=0 

e ^ ^ d n x , 

n=0 

6 ^ e n . 



-n-2 



(5.11a) 
(5.11b) 
(5.11c) 



n=0 



Note that Co = 1. Asymptotic solutions to Eqs. (|5.8p can be made by substituting the above 
expansions into these equations. It turns out that the following pairs of sums result in solutions to 
Eqs. (|5.8p in the limit x — > oo: 



oft j 0*0 



m=0 
„2x 



Cui—m 



X 



Cui—m—2 



(5.12a) 
(5.12b) 



m=0 
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and 

= z~ 2x E <^" C "" m " 2 , (5.13a) 

m=0 

ftjto = Y. h ™ x ~ c "~ m - ( 5 - 13b ) 

m=0 

It is easily shown that Oq and 6q are both arbitrary constants whose values determine the values 
of the remaining coefficients. Note that we have chosen the solutions for a q u ■ and ■ so that the 
mode functions F q ■ and ■ are well behaved in the limit r — > oo, which is the correct boundary 
condition for these modes. 

The boundary conditions for most of the pairs of sums are now easily obtained. First note that 

lim a p , Ax)a q , Ax) = alal lim x~ 2 e' 2x = , (5.14a) 
Jim c? Uti {x)lP Uti {x) = ogftg, (5.14b) 
lim al Ax) Pi Ax) = lim x" 4 = . (5.14c) 



The Wronskian condition (|5.7p then gives 

lim a p ujJ (x)^ J (x)=a p bl = l. (5.14d) 

x— >oo 1J u 

Note that all that is uniquely fixed is the product of the two functions in this limit. However, this 
is all that is necessary for the computation of {T a ^). Solving Eq. (|5.10b|) using (|5.14dp one finds 
that to first order in Cuj that 

<i(^)^ >i (^) = l (5-15) 

for all values of x. 

To obtain a final boundary condition, we must impose the constraint that the p-solution be 
well-behaved at the inner boundary. This is difficult because we do not in general have knowledge 
of the metric in this region. We anticipate that in the region r ~ max(|A|, |B|), the metric will vary 
significantly from its large r approximation (|2.2p . Let us choose a distance a which is sufficiently 
larger than max(|A|, |B|) so that the approximation (|2.2p can be assumed accurate for r > a, but 
which is small enough that a r e , with r e the value of r at which we wish to evaluate the stress- 
energy tensor. Because we remain agnostic about the metric for r < a, we cannot know the relative 
size of the K and I contribution to the p- modes at r = a. Looking at Eq. (|5.6bp . one might expect 
that generically 

"^■(^oKj+i^a) ~ Pl,j{ua)K j+1 {Lja) . (5.16) 
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We anticipate that the most important contributions to the stress-energy tensor at the radius r e 
come from to ~ r" 1 . Since a <C r e , we therefore expect that wa will be small. At small x, Kj{x) 
and Ij(x) behave as x~i and x J respectively. Using Eq. (|5.16p . it follows that 

^ 3 (ua)~{ua?i+ 2 ai.{ua). (5.17) 

Multiplying both sides by ■ (wa) and using Eq. (|5.15p , we then have 

^HflSjM ~ M 2 i+X>a)/^M ~ ^ • (5-18) 

Note that the smallest j we are interested in is j = | , and hence Eq. (I5.18P falls at least as fast as 
r~ 3 . Recall that we are only interested in keeping terms to order toC ~ C/r e . For sufficiently large 
r e , Eq. (|5,18p will always be negligible compared to C/r e , and hence we can treat this as zero, so 
that 7 

lim^ J (x)^ i .(x) = 0. (5.19) 

Solving Eqs. (|5.10p subject to the boundary conditions Eqs. (|5.14ap . (|5.14cj) . (|5.14d[) . and (|5.19p 
is now straightforward: 

POO 

= Coj / [K](y) - K? +l (y)] dy , (5.20a) 

J X 

<M0LAx) = 1, (5.20b) 



= 0, (5.20c) 



/^■M,^) = y - ^+i(y)] • (5.20d) 

To compute the quantities A\ and A2, first solve Eqs. (|5.2[) for F w j and G u j in terms of F U j 
and G^j and then substitute into Eqs. (|4.14p . Next use Eqs. ()5.6p to find expressions for A\ and 
A2 in terms of ctj and /3j. Then use the explicit forms for the products in Eqs. (15.20p . Expanding 



It is conceivable that for certain pathological spacetime metrics this relationship will not hold, and hence the 0? 
will be important. That is why our argument is not completely general. 
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to linear order in C we find 

oo 



^ r ( roo 

+ Cx [if 2 (x) - tf? +1 (s)] f [/ 2 (y) - J? +1 (y)] dy - - [^(x^x) + 7 i+1 (x)l^ +1 (x)] 

Jo w 

+ B(2j + 1) [I J (x)K J (x) - I j+1 (x)K j+1 (x)}\ + ^ + A 

00 r C /"oo 

A 2 = ^2x(j+±) 2 |c^(x)I i+1 (x) ^ [K 2 (y)-K 2 +1 (y)]dy 



■** 1 A 



(5.21a) 



Ci^(x)K J+1 (x) / [lj(y) - if +1 (y)] dy - -I J+1 (x)^( 



+ B [I j (x)K J (x) - I 3+1 {x)K j+1 {x)] } + — + Ax 



+ — + iAx . 
4a; 1 



(5.21b) 



The sums over j for each term in A\ and ^2 have been computed in the Appendix. The integral 
terms in the formula for A\ are given by Eq. (|A.5ip . while the other terms in the sum are given by 
Eqs. (|A.17|) and (|A.22p . The integral terms in the formula for A2 are given by Eq. (TOOl) . while 
the other terms in the sum are given by Eqs. (|A.24|) and (|A.22p . The result is 



C I 2x 2 



A-2 



C[ fx 3 



dy 
V 

y 



e~ 2y + 2x - \ - xe~ 2x + \e~ 2x 



_ 2 . x 2 e~ 2x -J- ±™~ 2x -J- i^ _2:r ' 



£(j-x) + 2B(J-x) + | ; + jA f 



(5.22a) 



H - ^XG ~\~ ^G 



+ 2x 2 x 6 



-(k 2 -§x 3 ) 



+ 2£x(± - x) + — + ^Ax . 
Simplifying these two expressions, and replacing 00 with x/r, we find 



xr + 2Ax - Cxe~ 2x + \Ce~ 2x + 2Cx 2 



-2y 



(5.22b) 



(5.23a) 



A 2 = | x 2 r + 2Ax 2 _1 C _ 2 Cx 2 e -2 X + l Cxe -2x + l Ce -2x + 4^3 j°° ± e ~2y (f; 33^ 

J x y 

To complete the computation of {T a f,), first recall that the analytic contribution falls as r~ 6 at 
large r, and hence can be neglected. The numeric contribution, Eqs. (|4.ip . is thus the dominant 
contribution. Substituting Eqs. (|5.23p into Eqs. (14. ip and changing the integration variable from 
uj to x = wr gives to linear order in C/r 

C r °° 



{T\) 

(T r r) 

(T%) = (T%) 



7r 2 r 5 

C 
7r 2 r 5 
C 
2vr 2 r 5 



i 2; 2 e - 2 ^-x 3 e- 2a; + 2x 4 



00 ^e- 2 ^ 



tX e 



^ »Z< G dx j g 37 G 



3 x 



dx , 
"dy 



-2y 



2 3 -2a: _ 12 p -2x _ l_ p -2a: _ 4 4 



y 

dy 



-2y 



(5.24a) 
dx , (5.24b) 
dx . (5.24c) 
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Integrating the last term in each expression by parts we find that in the (f , r, 9, <p) basis, 



< T ^> = iTTTT dia g ( 4 > 2 ' ~ 3 ' " 3 ) • ( 5 - 25 ) 



VI. ANALYTIC COMPUTATION OF THE STRESS-ENERGY TENSOR FOR SCALAR 

FIELDS 

We now proceed to a computation of the quantities {(j) 2 ) and (T a b) f° r massless scalar fields. As 
can be seen from Eqs. (|2.4p . (|2.5p . and (|2.6p these quantities contain sums and integrals over the 
radial mode functions. Thus we begin by finding approximate solutions to the mode equation (|2.8p 
which are valid at large values of r. 

To solve the radial mode equation it is useful to first change variables in such a way that the 
modes and the equations they satisfy are similar to those which are solved for the spin ^ field. We 
begin by defining a new mode function F such that 

1 'f^F-B^). (6,a) 



r 1 / 2 \ 2r dr 
To first order in A and B the inverse relation is 

F^{ S+ B -^S + B f). (6,b, 

If ()6.1ap is substituted into the mode equation (|2.8p then the resulting equation has terms of zeroth, 
first, and second order in A and B. Since we are working to first order, the second order terms 
can be dropped. One of the first order terms is proportional to Bd 3 F/dr 3 . This term can be 
eliminated. To do so first take one derivative with respect to r of the zeroth order equation and 
solve it for d 3 F/dr s . The result is 

d 3 F _ 1 d 2 F | f 2 | (I + |) 2 ( 3 \dF (u? [ 5(t+j)* \ p 
dr 3 2r dr 2 \ r 2 2r 2 J dr I 2r 2r 3 J 

Substituting this into the equation for F which contains both zeroth and first order terms, one 
finds that to first order F satisfies the equation 



d 2 F 1 dF 
+ 



, (6.3) 



dr 2 r dr 

where, as before, C = A + B. 

There are two solutions to this equation, F^ e and F^ £ , for each value of ui and i. The Wronskian 
condition (|2.10j) becomes 

dF q dF p 1 

——F p - -——F q = -- . (6.4) 
dr dr r 
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Defining x = wr as before, and letting primes denote derivatives with respect to x, one finds that 
Eqs. (]6.3p and (]6.4p become 



F" + -F' 

x 



2Cu (l + h) 2 
1 + + - 



F = 0, 



(6.5) 



and 



pi'pp _ pp'pi 



1 

x 



(6.6) 



If C = 0, then the solutions to Eq. (|6.5p are again modified Bessel functions. As for the spin ^ 
field, we define new functions ay and /3j, where j = I + ^, so that 



F'(x) = a UJ j(x)l' j (x) + /3 UJ j(x)K' j (x) 



(6.7a) 
(6.7b) 



In general, this can always be done, since we can treat this as two equations in two unknowns. 
Demanding that these equations be compatible implies 



^,j(x)Ij(x) + ^ J (x)K j (x) = . 



(6.8) 



Substituting Eqs. (|6.7p into Eq. (|6,5p . and using the fact that the Bessel functions satisfy the 
modified Bessel's equation (jA.ip . we find 



a^ i (x)/j(x)+^ ii (x)^(x) 



2Cu 



n 



[x)Ij{x) + p u j(x)Kj(x)] 



(6.9) 



Solving Eqs. (|6.8p and (|6.9|) for a'^ ■ and and using the modified Bessel function Wronskian 
Eq. ()A.4p . one finds 

o/ UJ (x) = 2Coj[a U3jj (x)K j (x)Ij(x) + l3 u}tj (x)K?(x)] , (6.10a) 
p u j(x) = -2CLj[a U)j (x)l](x) + ^, j (x)K j (x)I j (x)] . (6.10b) 

The Wronskian (|6.6p in terms of the a's and /3's is identical to leading order to Eq. f|5.7j) : 

- <;(*)/^(z) = 1 • (6-11) 
As in Section [Vj we will work not with Eqs. (|6.10p . but with the derivatives of the corresponding 
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products. We find 



:x)al d (x)]' = 2Cu{2a? u>j (x)al tj (x)K j (x)I j (x) 

{ [<^)/^-(*) + /Sj^-W] 
+ 2/3^.(x)/3^.(x)^(x)/,(x)} . 



2Cw 
2Cw 
-2Cw 



(6.12a) 
(6.12b) 
(6.12c) 

(6.12d) 



As for the spin i field, we keep only terms to first order in Cuj in these equations. Thus we 
need to only keep the variables on the right side to zeroth order. As before, we expect that all of 
the products except for -(x)^ -(x) will vanish to leading order, and because of the Wronskian 
condition Eq. (|6.11]) . we expect this product to be one to leading order. Hence Eqs. f|6. 12j) can be 
simplified to 



a 



(6.13a) 
(6.13b) 
(6.13c) 
(6.13d) 



Asymptotic boundary conditions at infinity can be obtained by solving Eqs. (|6.10p in the large 
x limit. Using the usual asymptotic expansions for Bessel functions, we find solutions of the form 



w,3 
IP 



^ ] a m x 

m=0 



m=0 



X 



Coj— m— 1 



(6.14a) 
(6.14b) 



and 



a lA x "> 



-2.1 



E' 

m=0 



-Ccu—m—l 



-Cui—m 



(6.15a) 
(6.15b) 



m=0 



The g-solutions are chosen so that F q vanishes in the limit x — > oo. Using Eqs. (|6.14p and (|6.15p 
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along with the constraint (|6.1ip . we find the three boundary conditions: 

} i ^ a l,j( x ) a l,j( x ) = °i (6.16a) 
lim oP^Wfijix) = 1, (6.16b) 
lim al 4 {x)ff, 4 {x) = . (6.16c) 

Note in particular that the boundary condition (|6.16b[) implies that the solution to Eq. (|6.13bp for 
all values of x is 



a 



L-M,W = 1- (6-17) 



The fourth boundary condition requires special care in this case. It comes from the behavior of 
the p-solution, which must be regular at the inner boundary. As in Sec. V, we choose a distance 
a which is large enough that the approximation (|2.2p can be assumed accurate for r £ a, but 
which is small enough so that a <C r e , where r e is the value of r at which we want to evaluate 
the stress-energy tensor. We expect that at r = a the contribution of Ij and Kj to • will be 
comparable, so that 

a^a)/>a) ~ ^H^'W ■ (6-18) 
Since Ij(x) and Kj(x) for small x go as x 3 and x~ 3 respectively, this tells us that 

^.H-H^^H. (6.19) 

Multiplying both sides by /J'^wa) and using Eq. (|6.17p gives 

#>«)#jM ~ ^ ■ ( 6 - 2 °) 

We anticipate that the largest contribution to the stress-energy tensor at r = r e will come from 
oj ~ r^ 1 , and therefore expect that the right side of Eq. (|6.2U|) will be proportional to (a/r e ) 2 i . 
For j > |, this is negligible compared to Cu> ~ C/r e for large enough r e , and hence we treat it as 
zero. But for 3 = \ this term could be comparable to (or larger than) Cui. 
For j = 5 , it is useful to define the constant 7 such that 



7T 



^,1/2 



7 = lim 7" . (6.21) 

The factor of n is included to simplify later relationships. The limit au — > is appropriate, because 

we expect cj ~ r^ 1 , so this is equivalent to r e 3> a. Using Eq. (|6.17p we therefore find 

7LJ 7CJ 
lim /3 P , ,J(jja)B q -. ,„(wa) = — lim oP . ,Juja)B q -. ,A<joo) = — . (6.22) 
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This yields our fourth boundary condition 8 



(6.23) 



It is important to note that 7 is of linear order in the metric perturbations A and B. 

Using the boundary conditions (16. 16ft and (|6.23p one easily finds that the solutions to Eqs. (16. 13[) 



arc 



-2Cu; / Kj(y)dy 

J X 



70; 



<5,- 1/a - 2Coj \ If(y)dy . 
k Jo 



(6.24a) 
(6.24b) 
(6.24c) 
(6.24d) 



Starting with Eq. (|6.1aj) and using Eqs. (|6.7ap and (|6.24[) gives 



Pu:i(r)quje(r) 



1 A 



I j (x)K j (x) - 2Clo 



px poo 

K]{x) I]{y)dy + I](x) Kj(y)dy 

JO Jx 



+ 



J^ p -2x 

2x 



J, 1/2 



Boj 



(6.25) 



where the explicit form for Ky 2 , which is given in Eq. (jA.2h . has been used. 

To compute Si and 5s, one can first multiply Eq. (|6.25p by 21 + 1, subtract 1/r + A/r 2 , and 
then sum over £ with the result that 



E 

1=0 



(2£+l) Puje (r)q ul£ (r) 



1 A 

I2 



r 



+ ( i + 4 1 Prt-WJTjd) - 1] 



E 

i=i/2 
2Bjoj 



fX poo 

K){x) I I]{y)dy + I]{x) / Kj(y)dy 



[l j (x)K' j (x) + K j (x)l' j (x)]\+^-e- 2 > 



Using the sum formulas ()A. 13h and (lA.14j) along with the identity Eq. (lA.46j) yields 



E 

e=o 



(2e+l)p w e(r)q u t(r) 



A 



r 
2C 



dy , , , -2.1- 



— ( 2x f -< 

lx y 

..,.2 r -< 2 " + - r n 



y 



x Ax Bbj 70; _ 2x , 
r r 2 r 2rx 



(6.26) 



(6.27) 



It is conceivable that for certain pathological spacetime metrics this relationship will not hold, and hence the /3J 
will be important for j > J. That is why our argument is not completely general. 
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Then Si and S5 can be found by substituting the above result into Eqs. (|2.5ap and (|2.5ej) . and 
replacing 00 with x/r: 



Si 

s 5 



1 



dx x 2 ( 2x 



47r 2 r 5 J 

4^3 / dx ( he -2 ■" 



_2C7x 2 r^L e - 2 y + Cxe 

Jx y 

2Cx 2 [°° ^-e-^ + Cxe~ 2 

Jx y 



-2x 



The remaining integrals are straightforward to compute. The results are 

3C + 5 7 



Si 

S 5 



1607r 2 r 5 
C + 3 7 



487r 2 r 3 

Next 5*4 can be found from £5 using the relation (I2.9ap with the result that 

C + 3 7 



S 4 



16vr 2 r 4 



(6.28a) 
(6.28b) 

(6.29a) 
(6.29b) 

(6.29c) 



The sum S3 can be computed in a similar way by first substituting Eq. (|6.25p into Eq. (|2.5c|) 
with the result that 



5, 



dw J 
4^1 ^ 



- + ^)\2j 3 Ux)K j {x)-AfCu 



Kj=l/2 
/•oo 

+ /J(x) / K 2 (y)dy 



Kj{x) \ I (y)dy 




,, 1 2Buj 3 rr , . 1 , , , , -1 urr 3Au/ 



70; _ 2x 2r 2 u 3 8Aruj 3 uj Auj Boo \ 
8rx 3 3 4 6r 3r J 

Next replacing oj with x/r and rearranging slightly one finds 



(6.30) 



S, 



1 



dx 



4 7r 2 r 3 ^ Q 

+ 2j 3 ^.(x).fc: j (x) + x] -Cx 



J2\(r + A) [2j 3 I 3 {x)K 3 {x) - j 2 + \x 2 ] - Bx[2j 3 I^K'^ 
■i=i/2 l 



4fKj(x) J I 2 (y)dy + 4j 3 I 2 (x) / K 2 {y)dy 



1 - :ir - frrx 3 - %Arx 3 + \rx + ±Ac - \Bx 



(6.31) 



Substituting Eqs. ([P9|> . fl30| . and (IQil) for the remaining sums and integrals yields 
1 



4 7r 2 r 3 



<M C 



( 3 x 2 x ) 



°° ^ p -2</ _ 2 3 2z 12 -2x _ ± Tp -2x 



(6.32) 



Finally, computing the remaining integrals using integration by parts where necessary gives 

15 7 - 7C 



S, 



960vr 2 r 3 



(6.33) 
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To find 5*2, it is easiest to use Eq. (|2.9b|h To leading order this relationship is 



2 2 dr r 1 r 2 4r 2 



(6.34) 



Substituting explicit forms from Eqs. (I6.29a|) . (|6.29bj> . (jfi.2Hcj) . and (IQ3l . we find 

9C + 257 



5 2 



1607r 2 r 5 



(6.35) 



It is now straightforward to substitute all of our expressions for S n into Eq. (12. 4j) to obtain the full 
expression for (T a b): 



(T a b ) 



1 



80?r 2 r 5 



{Cdiag (-1, 2, -3, -3) + 5 [eC + 3(£ - ±) 7 ] diag (2, -2, 3, 3) } 



(6.36) 



In general this result depends on the details of the geometry at small r through the parameter 7, 
but for a conformally invariant scalar field (£ = g) the 7-dependance cancels. We can also compute 
the quantity (4> 2 ) by using Eqs. (|2.6p and (|6,29b|h keeping in mind that the numerical contribution 
dominates the analytical contribution at large r, so that 

C + 3 7 



Note that unlike the stress-energy tensor, the value of (cp 2 ) depends on 7 for any value of £ 



(6.37) 



To evaluate Eq. (|6.36p or (|6.37p explicitly, we need to determine 7. To do so it is necessary to 
solve, analytically or numerically, the exact mode equation (j2.7j) when uj = i = for the p-function. 
This equation is 



d 2 Poo + /2 + J_ df_ _ _}_dh\ dm _ t RhpQQ = Q 
dr 2 \r 2/ dr 2h dr J dr 

In the asymptotic region, the metric becomes flat, and the solution takes the form 



Poo(r) = c + ci/r . 



(6.38) 



(6.39) 



Substitution into Eq. (|6,lb|) then gives the function F p , which we match to order 1/r with the 
general form Eq. (|6.7a,p : 



*S> = ^00(1 + ^ + 



ci B-A 



lim 

u>-H) 



cor 2r 



r l l 2 I lim 



<l/ 2 ( Wr ) [ — ) 



1 + ^1 • 

2r' 



(6.40) 
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Here the definition (|6.2ip has been used. Clearly 

/9/A 1/2 ~ 

= c , (6.41a) 



lim 



<l/ 2 (^ r ) { — ) 



2ci 

7 = — + B-A. (6.41b) 

Co 

Thus 7 can be found by solving Eq. (|6.38p either numerically or analytically and matching the 
solution with its asymptotic form (|6.39p to obtain the values of cq and c±. 

For example, consider the case of a Schwarzschild black hole, with metric functions given by 
Eqs. (|3.5p . so that A = B = M. Since the scalar curvature it! is zero everywhere, there are two 
linearly independent solutions to Eq. (|6.38p . one of which diverges at the horizon, and the other 
of which is constant everywhere. The latter is the p-solution, and we see from Eqs. (|6.39[) and 
(|6.41bp that 7 = 0. The resulting stress-energy tensor is given in (j3.8|) for all values of ^. Thus 
our use of the WKB approximation in Section IIIII gives the correct leading order behavior for the 
stress-energy tensor in the far field region. 

In contrast, consider a nonrelativistic star, by which we mean a static spherically symmetric 
weak field source whose stress-energy tensor is dominated by T l t = —p(r). Such a source will have 
R ~ 8irp. In flat space poo is equal to a constant. In the weak field limit we can therefore write 

Poo = c 2 + S P (6.42) 



and substitute it into the mode equation (J6J38J). Note that without loss of generality we can set 
5p = at the origin. Assuming that 5p is small and keeping only first order terms gives 

d ( 2 d5p\ 



i \ r-f = ^r 2 p{r)c 2 . (6.43) 
ar \ ar J 

For points outside the star this can be integrated twice to obtain 

2£Mc 2 

Sp = c 3 - , 6.44 

r 

where M is the total mass and | C3 1 <C \c 2 \- Matching to the form (|6.39[) , we see that cq = c 2 + C3 . 
Then to leading order, c\ = -2£Mc . Using Eq. (j6.41b|> with A = B = M gives 7 = -4£M. Thus 
to leading order 

^ 2 )sch = ^3 > ( 6 - 45a ) 
2 W r(l-6e), (6.45b) 



/ nr 



24 7r 2 r 3 



M 

(T a b ) Sch = {diag (-1, 2, -3, -3) + 5£diag (2, -2, 3, 3)} , (6.45c) 

(T a b )nr = {diag (-1,2, -3,-3)+ [5^- 30£(£- i)] diag (2, -2, 3, 3)} . (6.45d) 
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111 ], to leading order 



These are the results found previously in Refs. [<| and [hJ. As shown in 
at large r the quantity (cj) 2 } in Schwarzschild spacetime differs from its value outside of a static 
spherically symmetric star except in the case £ = 0, while (T a b) differs from its value outside of a 
static spherically symmetric star except in the cases f = I and f = 0. 

The argument that for £ = the asymptotic values of (cp 2 ) and (T" 1 ^) are the same outside a 
nonrelativistic star as for a Schwarzschild black hole can be generalized to the statement that for 
£ = the asymptotic values of these quantities depend only on the leading order geometry in the 
far field region for any star or black hole. The exact solution is in this case simply 

PooO) = c , (6.46a) 
1 f°° dr' lh(r') . A , , 

Here it is clear that poo must be the constant solution because it is the only solution that does not 
diverge at r = in the case of a star, or at the event horizon in the case of a black hole. Then 
matching to Eq. (|6.39p gives c\ = 0, and hence 7 = B — A. 

The above result, while quite general, does not hold for a wormhole. By a wormhole we mean a 
metric such that at some throat radius b, h(r) has a simple pole while f(r) is finite and well-behaved. 
This represents a connection between two asymptotically flat universes, which for simplicity we will 
assume have identical metric functions h(r) and f(r). The radius r = b is the minimum radius for 
which the metric is defined and represents the point at which the two universes join. For minimal 
coupling, £ = 0, the p solution will be the one that vanishes at infinity in the other universe, and 
therefore it is not the constant solution of Eq. (|6.46ap . but instead the solution that behaves like 
r _1 at large r in that universe. Note that unlike the case when there is an event horizon, J y/h/ fdr 
is finite at the throat of the wormhole. The exact solution to Eq. (|6.38p will satisfy 



dpoo c h(r) 

^ = ± -^W)' ( } 

where c is an unknown constant. Care must be taken at the coordinate singularity r = b. A 
temporary change of coordinates to r = b + s 2 shows that the singularity in Eq. (|6.47p is integrable, 
so that poo is continuous, but dpoo/dr changes sign as you cross the throat. In Eq. (|6.47|) . the ± 
refers to which universe you are in, with the + sign corresponding to our universe and the — sign 
to the other universe. Keeping in mind that poo(°°) = in the other universe, one can integrate 
from 00 in that universe to the throat b of the wormhole and then in our universe integrate out 
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from the the throat to a radius r P with the result that 



, . f b dr h{r) r°dr h{r) „ f°° dr h{r) f°° dr h{r) fB ao . 

If r e is chosen so that it is large enough that we can treat the metric as flat for r > r e then 



f°° dr h(r) c 
P00 (r e ) = 2c J^l . 6.49 

Comparison with Eqs. (16. 39ft and (|6.41b|) shows that, for a scalar field with minimal curvature 
coupling £ = 0, 

-i 



7 



dr I h(r) 



+ B- A. (6.50) 



VII. NUMERICAL COMPUTATIONS 

There are several cases in which we have done numerical computations of the stress-energy 
tensor for massless fields far from the event horizon of a black hole or the throat of a wormhole. 
These are important because they are computations of the full renormalized stress-energy tensor 
for the quantum fields and thus serve as important checks on the WKB calculation and the exact 
analytic calculations in the previous sections. 

The first numerical computation of the stress-energy tensor for a scalar field in the Boulware 
state was done in the case of conformal coupling in [15(. However, it is not possible to deduce 
the large r behavior of the stress-energy tensor from the plots in that paper. In [1C| some results 
of numerical calculations for scalar fields with arbitrary coupling £ to the scalar curvature were 
shown. As discussed in that paper, the results agree with those in Eq. (|6.45cj) to within at least 
two digit accuracy. 

For the massless spin 1/2 field we have computed {T a b) in Schwarzschild and extreme Reissner- 
Nordstrom spacetimes. It is clear from the plots in Fig. 1 that the large r behavior of the compo- 
nents goes like M/r 5 . Since the metric for the extreme Reissner-Nordstrom spacetime is 

+ ^ (7.1) 

the values of A and B in Eq. A2.2H are A = B = M, the same as for Schwarzschild spacetime. Thus 
the leading order behavior of the stress-energy tensor predicted by Eq. (I5.25h for both Schwarzschild 
and Reissner-Nordstrom spacetimes is 

(n> = 40^ diag(4 ' 2 '" 3 '" 3) • (7 - 2) 
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FIG. 1: Numerical calculation of the stress-energy tensor for the massless spin 1/2 field in 
Schwarzschild (left) spacetime and the Extreme Reissner-Nordstrom (right) spacetime with the 
renormalization parameter /j, = M _1 in the notation of Ref. In each plot, the solid curve is 
(r*t), the dotted curve is (T r r ), and the dashed curve is (T e g) = (T^). 



/ 


h- 1 


A 


B 


3207r 2 r 5 (T a 6 ) 


1 


1 - b/r 





6/2 


26diag(4,2,-3,-3) 


1 - 6/2r 


1 - b/r 


6/4 


6/2 


36diag(4,2,-3,-3) 


l-b/r + b 2 /2r 2 


1 - b/r 


6/2 


6/2 


46diag(4,2,-3,-3) 



TABLE I: Leading order stress-energy tensors for three wormhole metrics. 



Fitting the numerical data to a series in inverse powers of r, we find that in each case there is 
agreement with the results in Eq. (j7.2|) to at least two digits. 

Numerical computations have also been done for three different wormhole metrics. In each case, 
the throat is at r = b. The explicit form of the metric functions / and /i~ 1 , together with the 
corresponding coefficients A and B and the predictions of Eq. ()5.25j) are given in Table 1. The 
stress-energy tensors for these metrics are shown in Fig. 2. It is clear from the figures that the 
large r behavior of the stress-energy again goes like 1/r 5 . Fitting the numerical data to a power 
series in inverse powers of r we find in each case agreement with the analytic results displayed in 
Table 1 to within approximately two digits. 
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FIG. 2: Numerical calculation of the stress-energy tensor for the massless spin 1/2 field in three 
wormhole metrics with renormalization parameter [i = b 1 in the notation of Ref . [14| . All three 
have /i(r) _1 = 1 — b/r. The metric function f(r) is given by /(r) = 1 (top left), f{r) = 1 — b/2r 
(top right) and f(r) = 1 — b/r + b/2r 2 (bottom). In each plot the solid curve is (T* t ), the dotted 
curve is (T r r ), and the dashed curve is (T e g) = (T"^). 



VIII. COMPARISON WITH PREVIOUS CALCULATIONS FOR SCHWARZSCHILD 

SPACETIME 

It is possible to compare our results for both the scalar and spin 1/2 fields with previous 
computations of the leading order corrections to the gravitational potential <3? in Schwarzschild 
spacetime which have been found by computing one loop corrections to the graviton propagator. 
Such corrections have been computed by Duff [l] for both the massless conformally coupled scalar 
field and the massless spin 1/2 field and by Hamber and Liu [jj] for the massless minimally coupled 
scalar field. 
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First note that the one loop correction to the graviton propagator for a given type of quantum 
field gives the leading order correction to the gravitational potential <I> from that field [J. Here <3? 
is defined by the relation 



9tt 



(1 - 2$) . 



1.1) 



Clearly at the classical order $ = M/r. 

Another way to obtain quantum corrections to $ is to solve the linearized semiclassical backre- 
action equations. This has been done previously for massless scalar fields with arbitrary coupling to 
the scalar curvature in Schwarzschild spacetime [l(J and in the exterior region of a nonrelativistic 
static spherically symmetric star [J]. If the metric functions in Eq. (|2.ip are written as 

2m(r) 



f(r) 
h(r) 



r 

2m(r) 



e 
n -l 



then the semiclassical backreaction equations take the simple form 

dm(r) 



dr 
#(r) 
dr 



4irr 



47rr 2 (r* t ) , 

(T\) - (T't) 
1 — 2m(r) jr 



3.2a) 
3.2b) 

3.3a) 
3.3b) 



Substituting Eq. (|6.45cf) into Eqs. ()8.3|) and integrating gives 



mr 



<Kr) 



M 



1 



;i - ioo 



207rr 2 
M(3 - 20£) 



which in turn implies that 



M 



+ 



1 

45-7T 



307rr 3 



1\ 1 



6 / 67T 

For £ = 1/6 this gives the same value as that found by Duff 111 and for f = it gives 

n □ 

found by Hamber and Liu [3J. The same agreement was found previously in Ref. [10] 
numerical results and the WKB calculation discussed in Section IIII1 

For the massless spin 1/2 field, substituting Eq. (|7,2p into Eqs. (|8.3j) and integrating gives 



(8.4a) 
(8.4b) 

(8.5) 

the value 
using our 



m(r) 



M 1 + 



1 



57rr z 



M 



3.6a) 
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and 



M 2M ( 1 . 



This value for the correction to <I> is twice as large as that found by Duff |l| for the spin ^ field. 
This is due to the fact that Duff considered a two component field and we are considering a four 
component one. 

One would expect the graviton propagator computations discussed above to give the correct 
changes in the gravitational potential in the weak field approximation, i.e. for the exterior of 
a static star. The reason that we get the same corrections for massless fields in Schwarzschild 
spacetime is that for the conformally coupled scalar field and the spin 1/2 field the leading order 
stress-energy depends only on the local geometry in the far field limit. As discussed previously, for 
the minimally coupled scalar field the leading order stress-energy at large r is the same for a static 
star as for a Schwarzschild black hole. 



IX. DISCUSSION 



In this paper analytic expressions have been computed for the leading order asymptotic behav- 
iors of the quantities ((ft 2 ) and (T a b) for arbitrarily coupled massless spin fields, and of (T a b) for 
the massless spin \ field, in the case that the fields are in the Boulware state. The results are 
expected to be valid for most static spherically symmetric spacetimes for which the metric has the 
asymptotic form (|2,2p . As discussed in detail in Sections IVl and IVTl there may be some spacetimes, 
which are probably pathological, for which our results are not valid. For this reason it is important 
to compare these analytic results with those obtained in other completely independent calculations. 

For massless scalar fields with arbitrary coupling £ to the scalar curvature in Schwarzschild 
spacetime a combination of the WKB approximation for the modes and conservation of the stress- 
energy tensor has been used in Section III to compute the asymptotic behaviors of the quantities 
(</> 2 ) and (T a b ). The results are in agreement with the more general analytic calculations mentioned 
above. 

The asymptotic behaviors of {(j) 2 ) and (T a b) in the region outside of a nonrelativistic, static 
spherically symmetric star have been computed analytically using different methods than ours in 
Ref. [4]. Their results arc in complete agreement with ours. 

The asymptotic behavior of the gravitational potential as defined in Eq. (|8.ip has been obtained 
by computing one loop corrections to the graviton propagator for massless scalar fields with con- 
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formal [l| and minimal coupling [3J] as well as the massless spin | field Another way to obtain 
the corrections to the gravitational potential is to solve the linearized semiclassical backreaction 
equations using the analytic results for the asymptotic behavior of the stress-energy tensor. This 
has been done previously for massless scalar fields with arbitrary coupling to the scalar curvature 
in Schwarzschild spacetime 10] and in the exterior region of a nonrelativistic static spherically 
symmetric star 4\. In both cases the results agree with the one loop calculations in [1, 3]. This is 
to be expected since there is no difference in the leading order behavior of the stress-energy tensor 
in Schwarzschild spacetime versus that of the exterior region of a nonrelativistic static spheri- 
cally symmetric star for the cases of conformal or minimal coupling In Section VIII we have 



solved the linearized semiclassical backreaction equations for the massless spin i field and found 

[ 

agreement with the corresponding one loop calculation in [1|. 

Numerical computations of both ((f) 2 ) and (T a b) have been made for massless scalar fields with 



In these calculations the 



arbitrary coupling to the scalar curvature in Schwarzschild spacetime 
full stress-energy tensor has been computed. This has also been done here for the massless spin ^ 
field in Schwarzschild spacetime, the extreme Reissner-Nordstrom spacetime, and three wormhole 
spacetimes. The results are displayed and discussed in Section VII. By fitting the numerical data 
to power series in 1/r we have found agreement with the analytic results discussed above. 

As mentioned in the introduction, the stress-energy tensor for quantized fields is intrinsically 
nonlocal. An interesting question that was raised in the comparison between static stars and 



Schwarzschild black holes in Ref. 



Ill ] is whether the leading order behavior of the stress-energy 



tensor in the asymptotic region of a static, spherically symmetric, asymptotically flat spacetime is 
local or nonlocal. As stated in the introduction, by local we mean that the leading order terms in 
the stress-energy tensor depend only on the leading order deviations of the metric from a flat space 



metric. In Ref. [ill ] it was found that when the geometry at large values of r is the Schwarzschild 
geometry, the leading order behavior of the stress-energy tensor for massless scalar fields is nonlocal 
except in the special cases of conformal and minimal coupling. This result was extended in Ref. Q] 
to the D > 4 dimensional Schwarzschild- Tangherlini geometry. 

Examination of Eqs. (I5.25P and (16.36P indicate that in the asymptotic region the stress-energy 
tensor is always local for the massless spin ^ field and for the conformally coupled massless scalar 
field. We believe that these results combine with those of Ref. [5] to give strong evidence that the 
locality of the asymptotic behavior of the stress-energy tensor is a consequence of the conformal 
symmetry exhibited by these fields. For the case of the minimally coupled massless scalar field 
the leading order asymptotic behavior of the stress-energy tensor is local for any static spherically 
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symmetric spacetime which contains either a star or a black hole, but it is nonlocal if the spacetime 
contains a wormhole. Thus in many cases of physical interest the leading order behavior of the 
stress-energy tensor in the asymptotically flat region depends only on the local geometry there. 
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In Sections [V] and IVII expressions for the stress-energy tensors of massless spin 1 /2 and spin 
fields are derived which contain sums of products of the modified Bessel functions Ij(x) and 
Kj(x), with half-integer values for the index j. In this Appendix we derive identities which make 
it possible to compute those sums. Throughout it is assumed that j takes on only half-integer 
values. The functions Ij(x) and Kj(x) are linearly independent solutions to the equation 
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Appendix: Bessel Function Identities 




(A.l) 



For j 




they are 



1-1/2 ( x ) 



2 cosh x 



\/2ttx 




2 sinh x 





(A.2) 



All other half-integer values can be computed using the recursion relations 



K j±1 (x) = -K' j {x)± 3 -K j {x) . 



(A.3) 



These functions also satisfy the Wronskian condition 



I j {x)K' j (x) -Kj{x)l'j{x) 



1 



(A.4) 



x 
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which, with the help of the recursion relations, can be rewritten as 



x [I j {x)K j+1 (x) + Kj{x)I j+ i(x)} = 1 



(A.5) 



For < x <C v7 an d j > 1, 



l s {x) 



1 



x\i 



r(j+i) V2 



2\ 1 



while for i>j' 



2 V x 



1 



K,(x) 

One of the sums which we need is 



\/2irx 

7T 



1 + o ( — 
1 + OI j 



1 + 



(A.6) 



1 + 



(A.7) 



Q x = lim [I i (x) J ft: i _i(x) + I i+ i(x)K J -(x)] 

i=l/2 



(A.8) 



To evaluate this sum first multiply by x, then take the derivative with respect to x and use the 
recursion relations (|A.3|) to obtain 



— ( 



Qi) = lim V] x [I 7 -_ 1 (a;)ifj_i(x) - Ij +1 (x)K j+1 (x)] 

J— >oo ' — » 



(A.9) 



i=i/2 

Most of the terms cancel, as can be seen by shifting the indices on the two sums. The result is 

J-i J+i 



-f-(xQi) = lim 

ax j^oo 



xIj(x)Kj(x) — xIj(x)Kj(a 

j=-l/2 j=3/2 



xI_ 1/2 (x)K_ 1/2 (x) + xI l/2 {x)K 1/2 (x) - liir^ [xlj + i(x)iirj + i(x) + x/j(x)iirj(x)] 



(A.10) 



Using Eqs. (|A.2j) and ()A.6[) one finds that 

" 1 • 

Integrating this equation from to x, and dividing the result by x yields the identity 

oo 

Qi=J2 [I j (x)K j - 1 (x)+I j+1 (x)K j (x)] = l. 
i=i/2 



(A.ll) 



(A.12) 
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Using the recursion relations (|A.3P and the Wronskian condition (1A,4|) then gives the closely related 
identity 

oo 

Y [2jI j (x)K j (x)-l] = -x. (A.13) 
3=1/2 

A third identity can be found by taking the derivative to obtain 

oo 

Y 2j[P j (x)K j (x) + Ij^Kfa)] = -1 . (A.14) 
3=1/2 

Another sum that we need is 

J 

Q 2 = lim iU+h) [Ij^Kjix) + I 3+ i{x)K ]+l {x)\ - l} . (A.15) 

3=1/2 

It can be obtained by first shifting the indices on the second term and then combining the two 
terms so that 



Q 2 = lim < 

J— >oo 



J J+l 
Y [U+l)Ij(x)K j (x)-l]+ ^O'-iKjf^W 

U=l/2 J= 3 / 2 

lim { V [2jI J (x)K 1 (x)-l} + (J+^Ij +1 (x)K J+1 (x)\ . (A.16) 

J— >oo I * — ' I 

3=1/2 



The limit of the sum is given by Eq. (|A.13|) . and the limit of the last term can be evaluated using 
Eqs. (|A.6p . with the result that 

oo 

Q2=Y (O'+l) [Ij(x)K,(x)+I J+1 (x)K J+1 (x)} - 1} = i - x . (A.17) 

3=1/2 

A slightly more difficult sum is 

J 

Q 3 = lim Y {(J+l) 2 ViWKji?) ~ I 3+ i{x)K 3+1 {x)} - ±} . (A.18) 

J— >oo * — * 

3=1/2 

As before we begin by shifting indices and combining terms, so that 
J 



Q 3 = lim { Y [(j+ 1 2) 2 I^)K 3 (x)-(j-l) 2 I J (x)K J (x)- §] - (J+\?Ij +1 {x)K J+l {a 

U=l/2 

(A.19) 

The limit of the last term must be handled carefully. With the help of the recursion relations (|A.3j) . 



36 

the Wronskian ()A.5|h and the limiting forms (|A.6p . this term can be rewritten as 

{J+\) 2 I J+1 {x)K J+l (x) = X !f { j^ llj(x) - Ij+2{x)\ K J+1 (x) 

= \ (j + [1 " xI J+1 (x)Kj(x) - xI J+2 (x)K J+l (x)} 

J 

= i j + o ( j-i) °( J_1 ) • ( A - 2 °) 

3=1/2 

Substituting this into Eq. (1A.19P gives 

J 

Q 3 = hm V [2jI J (x)K j (x) - 1] + ± . (A.21) 

J— >oo ' — * 

3=1/2 

Then using the identity ([A7T3T) one finds that 

oo 

Q 3 = E {^'+l) 2 ft^)^^) " Ij+i(x)K j+ x(x)] ~\} = \-x. (A.22) 
3=1/2 

Multiplying this sum by x and using the recursion relations (|A.3P to rewrite /j in terms of Ij+i, 
and i^j+i in terms of K 3 yields 

oo 

£ {(i+i) 2 [x/j +1 (x)i^(x) + x/ j+1 (x)i^(x) + I 3+1 {x)K 3 {x)] - \x) = \x - x 2 . (A.23) 

3=1/2 

The left side is a total derivative, so we can integrate and multiply by a factor of 2 to obtain the 
identity 

oo 

£ {2{j+\fxI j+1 {x)K 3 (x) - ±x 2 } = \x 2 - fx 3 . (A.24) 

3 = 1/2 

Another quantity we need is 

oo 

= £ [j 3 /,(x)i^(x)-j 2 + ±x 2 ] . (A.25) 

3=1/2 

Its value can be derived from Eq. (|A.24|) using Eq. (|A.5P in a slightly different way: 

oo 

£ {(j 2 + 2j + \)xI j+1 {x)K 3 (x) + (j 2 - |) [1 - xI 3 ix)K j+1 (x)] - ±x 2 } = ±x 2 - fx 3 . (A.26) 

3=1/2 

Writing the infinite sum as the limit of a finite sum and then shifting the index so that j — >■ j — 1 
on the first term yields 
J 

E {if - i) K(X)^_!(X) + 1 - Xl,(x)i^ + 1 (x)] - |X 2 } 



lim 

J— >oo 



■3=1/2 

+(J 2 + 2J + f)x/j +1 (x)iO(x) 



±x 2 -§x 3 . (A.27) 
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Next use the recursion relations ()A,3P to eliminate Kj±\ and use the limiting form (|A.6|) on the 
last term to obtain the identity 

oo 

£ {{j 2 - \) [1 - 2jI j (x)K j (x)] - \x 2 } = -fx 3 . (A.28) 
3=1/2 

Using the result (|A.13P gives the identity 

oo 

Q4 = [iflj^Kjix) - j 2 + \x 2 } = lx 3 - \x . (A.29) 

3=1/2 

Another identity is obtained from the derivative of this equation: 

[ZfljWKjix) + 2fl j (x)K' J {x) + x]=2x 2 -\. (A.30) 

3=1/2 

It is also necessary to evaluate integrals of products of Bessel functions and in some cases both 
integrals and sums of products of Bessel functions. One such integral is 

poo 

Q 5 = (2j+l) / [Kj(y) + K] +1 (y)] dy . (A.31) 

J X 

Using the recursion relations (|A.3P to rewrite 2jKj(y) in terms of Kj(y) and Kj + \{y) and (2j + 
2)Kj + i(y) in terms of K'- +l {y) and Kj(y) one finds 

/•oo 

Q 5 = / [Kf(y) + 2yK j (y)K' j ( y ) - K 2 +l (y) - 2yK j+1 (y)K' j+1 (y)] . (A.32) 

J X 

The right side is a total derivative, so 

/•oo 

Q 5 = (2j+l) / [K 2 (y) + K 2 +1 (y)] dy = xK 2 +1 (x) - xK 2 (x) . (A.33) 

J X 

The identity 

(2j+l) f [I 2 (y) + I 2 +1 (y)} dy = xl 2 (x) - xl 2 +1 (x) - -^_ 1/2 , (A.34) 
Jo 7r 

which is valid for j > — can be derived in a similar way. The last term comes from the lower 
limit in the special case j = —\. 
Next consider 

®® /*oo px 

Q e = J2 4 i J j( x ) / K 2 (y)dy + K 2 (x) I 2 (y)dy . (A.35) 

,=1/2 L Jo J 

To evaluate this quantity first take the product of Eqs. ()A.34p and ()A.33p and sum over j to obtain 
the identity 



J ( poo 

Jim * - J?+iO«0] / + K !+i(y)} dy 

+x [K 2 (x) - K 2 +1 (x)] j* [I 2 (y) + I 2 +1 {y)] dy} = 0. (A.36) 
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Then make the index shift j — )■ j — 1 on each of the second terms in the integrals. This adds 
additional terms at the lower and upper limits of the sum with the result that 
J 

■ I— : 



J ( poo px 

Urn \x[$-i(x)-*]+i(x)] J K 2 {y)dy + x[K 2 ^{x)-K 2 +l {x)] J I 2 (y)dy 

3 = 1/2 



r -, poo r -. px 

= x\l 2 _ l/2 {x) - I 2 l/2 {x)\ j K 2 l/2 {y)dy + x\K 2 _ l/2 {x)-K 2 l/2 {x)\ J I 2 l/2 {y)dy 

f poo px ~\ 

+ x\im\[lj +1 (x)-l}(x)] J Kj +1 (y)dy+[Kj +1 (x)-Kj(x)] j lj +1 (y)dyj .(A.37) 

Because we are taking the limit J — > oo, we can use the approximate Eqs. (|A.6|) to evaluate the I 2 
integral: 



ij(y)dy 



dy /y\ 2J 
r 2 (j+i) \2J 



1 + 



2 /x\2J+i 
(2J+l)r 2 (J+l) \2J 



1 + 



(A.38) 



Eq. ()A.38P is valid not only for large J, but whenever x <C \fj . 

A similar identity can be derived for the K 2 integral, but because the integral extends to infinity, 
we must only use the approximate Eqs. (|A.6[) up to some limit A, where A <C \[~J . Hence we first 
rewrite the integral as 



Kj(y)dy 



rV) 



r 2 (j) (2 



dy - 

\y 

2J-1 



13 



1 + 0| T 



+ 



l + O [ — )+() 



A 
X \ 

a) 



Kj(y)dy 



2J-1 



4J-2 \x 

Then the last term can be bounded using Eq. fA33l) with j = J - 1, so that 

f'QO A /*O0 

J K 2 (y)dy = —[KSW-Kj^A)]- J Kj^dy 



K 2 (y)dy . (A.39) 



< 



A 



T *j(A) 



2J- 



r 2 (j) /2 



4J-2 Vx 



2J-1 



a; 



2J-1 



A 2 

, + „i T 



(A.40) 



If we choose A noticeably larger than x, but still small compared to yj, then we can neglect both 
contributions that are suppressed by 0[(x/A) 2J ~ 1 ], so we simplify Eq. (EA391) to 



Kj(y)dy 



r 2 (j) 



i + o (y 



(A.41) 



(4J-2)(x/2) 2J - 1 

Again, Eq. ()A.41|) is valid not only for large J, but for any x <C \/~J if J > |. Combining Eqs. (|A.38j) 
and (lA~4ll with Eqs. (TATH]) gives 

Jim | - I 2 {x)] K 2 J+l {y)dy} = hm | [i^+i(*) - i^(z)] jf = . 

(A.42) 
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Substituting these two limits into Eq. (|A.37|) yields 



J ( poo px 

lim [ x VU^) - I hi( x )} J K !(y) d y + x i K U( x ) - K !+i( x )] J o 



l]{y)dy 



3=1/2 



K l/ 2 (y) d y + x 



K 2 _ 1/2 (x) - K 2 /2 (x) 



I 2 1/2 (y)dy . (A.43) 



Using the recursion relations Eqs. (|A.3h to rewrite the factors on the left, and using Eqs. (|A,2h to 
rewrite the factors on the right gives 



OO r n n rjp 

V 4 j IjWlfa) K 2 (y)dy + K j {x)K' j {x) I 2 {y)dy 

,=1/2 L Jx JO 

It is not hard to see that the left side is a total derivative: 
— ^ 2j I](x) / K 2 (y)dy + K 2 (x) / 

X j =1 / 2 L 



— < 



dy 



-22/ 



-2.y 



(A.44) 



(A.45) 



Next multiply by two and integrate from to x. On the left side, it is not hard to show using 
Eqs. (|A.6|) . (|A.38|) . and (]A.41j) that in the limit x ->• all the terms vanish for j > §. The j = \ 
term can also be shown to vanish in this limit, so we find 



I 2 Ax) / K 2 {y)dy + K 2 {x) I 2 (y)dy 



2x 



Jx y 



More complicated expressions can be evaluated using Eq. (|A.46p . For example, 

J C poo 

& = J™, E 1 + C?+i> - I l+i( x )] / [ K ](y) - K !+i(y)] d y 

^°°i=l/2 l Jx 



+ (j+l)x[K 2 (x)-K 2 +1 ( X )] / [/|(y)-/| +1 (y)]dy 



(A.47) 



First shift the sums such that the factor in front of the integral is always I 2 {x) or K 2 (x). As usual, 
this will introduce terms coming from the shifted sum. The J = \ terms are proportional to j — \ 
and thus vanish. Eqs. COL COSjh and (TPT1) can then be used to simplify the terms at large J: 



Qi = lim 

J— >oo 



J ( poo 

£ 1 + x/?(x) / [2jK 2 (y) - (j+\)K 2 +1 {y) - (j-^K^y)] dy 

=1/2 ^ ' 

+ xK 2 (x) £ [2jl 2 (y) - (j+\)I 2 +l (y) ~ {j-h)l]-M] dy) 

- (J+±)x Uj +1 (x) J [K 2 {y) - Kj +1 (y)] dy + K 2 J+1 (x) J [l 2 (y) - lj +1 (y)] dy 

OO f poo 

£ |l + xlf(x) J [2jK 2 (y) - {j+\)K 2 +l {y) - (j-h)«J-M dy 

+ xK 2 (x) f [2jl 2 (y) - (j+\)I 2 +1 (y) - (j-l)lU(y)] dy)-\. (A.48) 



J=l/2 
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Using Eqs. (|X33l) and (|X34l) the expression can be rewritten as 

°° f POO PX 



i=i/2 



/J l + l] +1 (x)K](x) - l](x)K 2 +1 (x)] \ + ^K 2 /2 (x) - \ 



(A.49) 



The identity ()A.46P can be used to evaluate the terms with integrals. The rest of the terms in the 
sum can be simplified with the help of Eq. (|A.5p , Then using the explicit form for hi Eq. 
(IA31) gives 



C OO J l ^~ > 

Q 7 = 2x 2 / ^e- 2y + x- xe~ 2x + x ^ [^(a;)/^^) + 7 i+1 (x)i^(x)] + ±e" 2:r - \ . (A.50) 
V 3=1/2 

This can be simplified using Eq. (IA.12j) to yield 

00 r poo 

Qr = E 1 + ti + D x t 7 K x ) " / [ K » " K i+i(y)] ^ 

,=1/2 1 A 

+ (j+±)x [tf?(x) - i^ +1 (x)] jf [J?(y) - lf +1 (y)] dy 
= 2x 2 [°° ^e-^ + 2x - \ - xe~ 2x + \e~ 2x . 

J x y 

Another expression we will need to evaluate is 

00 r r poo 

Qs^Yl \^<j+\) 2 \lAx)I ]+ i(x) J [K 2 (y)-K 2 +1 (y)]dy 

- K 3 (x)K J+1 (x) £ [I 2 { y) - I 2 +l ( y) ] dyX + 



(A.51) 



i=i/2 



(A.52) 



Begin by taking the derivative of this expression. Then use the recursion relations Eq. (|A.3P to 
eliminate the derivatives of the Bessel functions. Next use Eq. (|A.5j) to obtain 

J J V f POO 

±Q* = A™ E 2 ^H) 2 \ * + / W(v) - tf+M d y 

-+°° i=1/2 L I 

+ x [K|(x) + jT [7|(y) - /J +1 (y)] + I j+1 (x)K j+1 (x) - Ij(x)Kj{x)\ + 1 . 

(A.53) 



Eq. (|A.22|) can be used to simplify the nonintegrated terms in the sum. Then shift the sums so 
the factor in front of the integral is always I 2 [x) or K 2 (x). This shift does not introduce any new 
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terms at small j, but it does at large J. The result is 

('). - lim } 

dx 



d f°° 

-Q 8 = jim Yl 2xI H x ) / {(j+h) 2 i K Hv) - K hi(y)) + H-hf [k]-M - K M } d v 

+ 2xJtJ(z) f {(j+i) 2 [Jj(y) - Jj +1 (y)] + {j-\f [I^M ~ #(v)] } <*V 
•/ii 

+ 2(J+I) 2 x|/ 2 +1 (x)^ [^(y)-^ +i(y) ] dy + jS: 2 +i(x) ^ [J2( y) _j2 +i(y) ] d2/ 



+ 2x 



(A.54) 



The large J terms, with the help of Eqs. (fOj) . (|X38]) and (|A4i~|) . can all be seen to vanish except 
for one, which can be computed with the additional help of Eq. (|A.34|) : 

2(J+±) 2 xK 2 +1 (x) f I 2 (y)dy = x{J+\)K 2 +l {x) \xlj(x) - xI 2 J+1 (x) - (2 J+l) [* % +1 (y)dy 
Jo L jo 

= { [i - x/ J+1 (x)^j(x)] 2 - x 2 k 2 +1 (x)/ 2 +1 (x) _ ( 2 j+i)xk3 +1 (x) ^ ij +1 (y)dy\ 

= J+i + 0(J- 1 ). (A.55) 



Here Eq. (|A.5j) was used to help approximate the leading term. Substituting this in Eq. (|A.54|) 

gives 

j C ^ r /-oo 

= E ^O*) y {(i+^)M^|(y)-^ + i(y)] + 0'-^) 2 [^i(y)-^I(y)]} d y 



ki=l/2 

+ 2x*T 2 (x) f {(j+i) 2 [/J(y) - J| +1 (y)] + (j-i) 2 [i?_ x (x) - / 2 (x)] } dy 
Jo 



+ 2x 



i 



+ (J+i)| 
(A.56) 



Next use Eqs. (|AT33l) and ([Oil) in such a way that only the terms I 2 {y) or K 2 (y) appear inside 
the integrals. The result is 

J 



d 
dx 



Qs = 2x — \ + lim 



J— >oo 



J ( poo px 



j=l/2 



+ x 2 (j+I) [/ 2 +1 (x)if 2 (x) -I 2 {x)K 2 +l {x)} +x 2 (j-\) [l 2 _ 1 {x)K 2 {x)-I 2 {x)K 2 _ l {x)\ 



(A.57) 



The remaining integrals are given by Eq. ()A.46j) . The two terms in the final line in the sum are 
equal and opposite save for a shift of index. Thus 

Q 8 = 4x 2 [°° ^ e - 2 y+Ax-2xe- 2x -\+ lim (J+±) {l + x 2 [/ 2 +1 (x)K 2 (x) - / 2 (x)^ 2 +1 (x)] } . 

Jx y J -*°° 

(A.58) 



d 

dx 
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With the help of Eq. (fO]l and then Eq. (EO|) . the remaining limit is 

lim (J+|) {1 + x 2 - } = lim (2 J+l)xI J+1 {x)Kj(x) = . 

J — too J~^oo 

(A.59) 

After integrating both sides of Eq. (jA.58j) from to x, one can use the limiting forms Eqs. ()A,6p . 
(|A.38p . and ()A,4ip to show that the left side vanishes at x = 0. The result is: 



00 r f r-oa 

Q 8 = Yl 2 <i+\? Uj(x)I j+ i(x) / [K](y)-K] +1 (y)]dy 
- Kj(x)K j+1 (x) -I 2 +1 (y)] dy} +x 

- - l r 3 /"°° ^y p -2y _ 2 2 -2! , l- p -2a; , 1 -2a , 9 2 _ 1 1 

g «X> I O n »aj C- | g iA> O I g C- | - ' ■ 2 g 

j x y 

The final quantity which we need to evaluate is 

OO r /*00 /*X 

Q 9 = V ±j 3 I?(x) K](y)dy + ifKf(x) l](y)dy- 



(A.60) 



i=i/2 



(A.61) 



To do so first rewrite Eq. (|A.60|) as the limit of a finite sum, and then shift j — )■ j — 1 on the terms 
with integrals that contain or to yield 



lim 

J— >oo 



J r poo 

2x/ J (x) [(j+±) 2 I i+1 (x) - (i-i) 2 /^!^)] / K, 

j=l/2 l 



f(y)dy 



+ 2x^(x)[(j-i) 2 K J _ 1 (x)-(j+i) 2 ^ +1 (x)] / / 2 (y)dy + x 



/* OO /* 3? 

2x(J+l) 2 Ij(x)I J+1 (x) K 2 J+l (y)dy + 2x(J+\fKj(x)K J+l (x) I 2 J+l {y)dy 

Jx JO 



4 -.3 
3 X 



dy 



-22/ _ 2 2 -2a; , l Tp -2x , 1 -2a , o 2 _ 1 1 



(A.62) 



The large J terms can be evaluated with the help of Eqs. (lA.6p . (IA.38p . and (|A.4ip . and the 
recursion relations ()A.3P can be used to rewrite this as 



°° r pod 

| [tjxIjWfyx) - (4j 3 + j)lf(x)] J K]{y)dy 

+ [AjxKj^K'jix) - (4j 3 +j)K 2 (x)} [ X l](y)dy 

Jo 



i=i/2 



+ X 



4 3 f°° ^y p ~2y _ 2 2 -2! , l_„-2a! , 1 -2a , o r 2 _ 1 
3 / y e 3 X e -r 3 xe 6 t 6 



(A.63) 



The terms with derivatives can be evaluated using Eq. (1X441) and the remaining terms proportional 
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to j can be evaluated using Eq. (|A.46|) with the result that 



i=i/2 



4j 3 /?(x) / K?(y)dy + IfK^x) I I){y)dy - x 



'l r _ 4 T 3\ f°° <^y p -2y , 2 2 -2a: _ l_ p -2a: , ± p -2x _ n 2 _ J_ 
2 X 3 X ) I „, r 3 J ' e 3 xe ^ 12 e ^ x 12 



(A.64) 
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